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The RICIS Concept 


The University of Ho us ton -Clear Lake established the Research Institute for 
Computing and information Systems (RICIS) in 1986 to encourage the NASA 
Johnson Space Center (JSC) and local industry to actively support research 
in the computing and information sciences. As part of this endeavor, UHCL 
proposed a partnership with JSC to jointly define and manage an integrated 
program of research in advanced data processing technology needed for JSC’s 
main missions, including administrative, engineering and science responsi- 
bilities. JSC agreed and entered into a continuing cooperative agreement 
with UHCL beginning in May 1986, to Jointly plan and execute such research 
through RICIS. Additionally, under Cooperative Agreement NCC 9-16, 
computing and educational facilities are shared by the two institutions to 
conduct the research. 

The UHCL/R1CIS mission is to conduct, coordinate, and disseminate research 
and professional level education in computing and information systems to 
serve the needs of the government, industry, community and academia. 
RICIS combines resources of UHCL and its gateway affiliates to research and 
develop materials, prototypes and publications on topics of mutual interest 
to its sponsors and researchers. Within UHCL, the mission is being 
implemented through interdisciplinary involvement of faculty and students 
from each of the four schools: Business and Public Administration, Educa- 
tion, Human Sciences and Humanities, and Natural and Applied Sciences. 
RICIS also collaborates with industry in a companion program. This program 
is focused on serving the research and advanced development needs of 
industry. 

Moreover, UHCL established relationships with other universities and re- 
search organizations, having common research interests, to provide addi- 
tional sources of expertise to conduct needed research. For example, UHCL 
has entered into a special partnership with Texas A&M University to help 
oversee RICIS research ant education programs, while other research 
organizations are involved via the “gateway" concept 

A major role of RICIS then is to find the best match of sponsors, researchers 
and research objectives to advance knowledge in the computing and Informa- 
tion sciences. RICIS, working jointly with its sponsors, advises on research 
needs, recommends principals for conducting the research, provides tech- 
nical and administrative support to coordinate the research and integrates 
technical results into the goals of UHCL, NASA/JSC and industry. 
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Abstract 

A robust Hoc control design methodology and its application to a Space Station at- 
titude and momentum control problem are presented. This new approach incorporates 
nonlinear multi- parameter variations in the state-space formulation of H <*, control the- 
ory. An application of this robust H x control synthesis technique to the Space Station 
control problem yields a remarkable result in stability robustness with respect to the 
moments-of-inertia variation of about 73% in one of the structured uncertainty direc- 
tions. The performance and stability of this new robust H controller for the Space 
Station are compared to those of other controllers designed using a standard linear- 
quadratic-regulator synthesis technique. 


‘Currently, Research Scientist at Dynacs Engineering Co., Inc., Clearwater, Florida. 
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1 Aerospace Engineer, Mission Planning and Analysis Division, Member AIAA. 
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1. Introduction 


The Space Station Freedom will employ control moment gyros (CMGs) as primary 
actuating devices during normal flight mode operation, and it will utilize the gravity- 
gradient torque for the CMG momentum management [1,2]. An attitude determination 
system of the Space Station will employ rate gyros and star trackers to compute the states 
of the vehicle for control purposes. Multivariable, periodic-disturbance accommodating 
controllers have been developed and are being considered for actual implementation to 
the Space Station Freedom [3-7]. . 

As illustrated in Fig. 1, the Space Station will be assembled and maintained using 
the Mobile Remote Manipulator System (MRMS) and its Mobile Transporter (MT). The 
MRMS/MT carrying a large payload will cause significant changes in the inertia property 
of the Space Station; consequently, it will affect the overall performance and stability of 
the control system. Study results on the effects of such MRMS/MT operations (e.g., a 
“bay” translation along the pitch axis and 180-deg slew maneuver about the pitch axis) 
can be found in [7]. The study results of [7] indicate that some form of adaptive or 
robust control with more than 50% inertia variation margins is necessary to account for 
the large changes in the inertia property caused by the motion of the MRMS/MT and 
its large payload. The study results also indicate that a high-bandwidth controller has 
unacceptable transient responses during the payload maneuvers. 

In this paper, a robust control synthesis technique based on Hqo control theory is 
developed and applied to the robust control design problem of the Space Station dis- 
cussed above. This new approach incorporates nonlinear multi- parameter variations in 
the state-space formulation of H ^ control theory [8-10]. An application of this robust 
Hoo control synthesis technique to the Space Station yields a remarkable result in sta- 
bility robustness with respect to the mo ments-of- inertia variation of about 73% in one 
of the structured uncertainty directions. Such a 73% inertia variation margin is rather 
significant compared to the margin of 44% of a typical linear-quadratic-regulator (LQR) 
design [3-7] with nearly the same control bandwidth as the robust H 0 Q controller. 

This paper is organized as follows. In Section 2 a robust full-state feedback control 
synthesis technique based on the Hoo control theory is presented, which exploits the 
concept of input-output decomposition of the uncertain system parameters [10-13]. A 
“full-state” feedback control is considered since the full states of the vehicle are available 
from an attitude determination system of the Space Station. A similar approach for the 
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dynamic compensator design can be found in [13]. The linearized equations of motion 
of the Space Station are reviewed in Section 3. A robust H <*> controller is synthesized 
in Section 4, with special emphasis on the input-output decomposition of nonlinear, 
uncertain multi-parameters of the system. 

2. Robust iToo Control Synthesis 


Background 

In recent years there has been a growing interest in robust stabilization and control 
based on control theory [8-10,13-15], and substantial contributions have already 
been made to the state-space characterization of the H control problems [8,9]. Most 
standard /foo-related control techniques are, however, concerned with the sensitivity 
minimization with respect to the external disturbances, and are not directly related to 
the structured parameter uncertainty. Recently, a new way of incorporating parameter 
uncertainty in the robust H ^ compensator design is developed in [10,13] by converting 
the parameter-insensitive control problem into a conventional Hoo problem. The state- 
space solution to a standard H control problem in [8,9] is then utilized by redefining 
the structured parameter variations in terms of a fictitious input and output. 

In this section, such a robust H control synthesis technique developed in [10,13] 
is reviewed with special emphasis on the new concept of “directional” parameterization 
of nonlinear, uncertain parameters. Only the “full-state” feedback control case is con- 
sidered here since the Space Station control problem does not need the consideration 
of state estimation. A more general case with dynamic compensation can be found in 
[10,13]. 

The H 0 o space consists of functions which are bounded and stable. The i^-norm of 
a real-rational matrix T(s) is defined as 

||T|U = sup{ ||T( S )|| : ReM > 0} 

= sup linn ii 

w 

- supff[T(ju;)] 

where <j[T{jus)\ denotes the largest singular value of T(ju>) for a given u. 


In this paper, we consider a linear, time-invariant multivariable system described by 


[8,9] 


i(t) = A x(t ) + B\ w(t) + B 2 u(t) 
z(t) = C x x(t) + D n w(t) + D u u(t ) 


(la) 

(lb) 


where se(t) is an n-dimensional state vector and is assumed to be directly measured, 
w(t) an mi-dimensional disturbance vector, u(t) an m 2 -dimensional control vector, and 
z{t) a pi-dimensional controlled output vector. 

The transfer function representation of this system is given by 


tn(s) 

u(s) 


( 2 ) 


*(*) = [ p n(s) P u (s) ] 
while the plant transfer matrix P(s ) is related to the matrices in Eqs. (1) by 

P(s) = C x (sl - A)" 1 [ B x B 2 ] + [ D n D n } (3) 


Internal Feedback Loop 

Consider an uncertain dynamical system described as 


x 

z 


A Bj B 2 
C\ D n D\ 2 


1 

X 


w 

- 

u 


(4) 


where C j, D n , and D X2 are not subject to parameter variations. The system matrices 
possessing uncertain parameters in Eq. (4) are linearly decomposed into an internal 
feedback loop [11,12,13] as follows: 


(5) 


where the first matrix in the right-hand side is the nominal system matrix and A e is the 
perturbation matrix defined as 



/ 




\ 

X 

X 

z 

= 

A 

B x 

D n 

b 2 

D\2 

+ A e l 

w 


\ 

. 


) 

u 


A A 

A e = 


A A ABi A B 2 
0 0 0 


( 6 ) 


Suppose that there are 1 independent parameters p x ,. • ■ ,pi and that their variations 
are bounded as pj < p, < p,-, or |Ap, | <1. If A e is Unearly dependent of each uncertain 
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parameter, then it can be decomposed as derived in [11,13]. However, A e may contain 
elements that are nonlinear combinations of Ap;’s. Variations in the uncertain matrix 
elements, Ae^, are represented in a functional form as 

Ae, = Ae,(Ap) for i = 1, • ■ • ,q 


where A p = [Ap 1? • • • , Ap;] T and q is the number of the uncertain matrix elements of 
A e . The perturbation matrix A e is then decomposed as 


A e 


M t 

0 


E [ N x N w N u ] = -MEN 


(7) 


where the columns of M and the rows of N span the columns and the rows of A e , 
respectively; and 


E = 

If Ae,’s are linear in Ap.’s, then 


Aei(Ap) 


( 8 ) 


0 Ae ? (Ap) 

he above input-output decomposition can be rearranged 
to become a rank-one input-output decomposition, which has been applied to a real- 
parameter variation problem [13]. 

Define the following new variables 


*p = 


N x 0 N w N u ] 


w v 

w 

u 


w p = - Ez v , 


(9a) 

(9b) 


then the perturbed system, Eq. (5), and the input-output decomposition, Eq. (7), can 
be combined as: 


X 


' A M x B x B 2 ' 


X 

z p 

z 

= 

N x 0 N w N u 
C\ 0 D n 


w p 

w 





u 


w p — — Ezp 


(10a) 

(10b) 


where w p and z v are considered as the fictitious input and output, respectively, caused 
by the plant perturbation; and E is considered as a fictitious, internal feedback loop 
gain matrix. 
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The above internal feedback loop representation of the plant parameter uncertainty 
becomes useful for stability /performance robustness analysis discussed later in this sec- 
tion. In fact, the parameter-insensitive control synthesis problem becomes a convensional 
H 0 o disturbance attenuation problem, which can be easily solved by using the state-space 
formulation of the Hoo control theory. 

Directional Parameter Variations 

A “hypercube” in the space of the plant parameters, centered at a nominal point, is 
often used as a stability robustness measure [16]. The robust control synthesis problem 
is then to find a controller which yields the largest hypercube that will fit within the 
existing, but unknown, region of stability in the plant’s parameter space. In [16], a 
computational method is developed, which exploits the mapping theorem and the “multi- 
linear” property of the plant’s uncertain parameters. However, as shown later in this 
paper, the Space Station has the uncertain moments of inertia which appear in the 
internal feedback loop gain E as nonlinear functions. 

One way to overcome the presence of such uncertain parameters in the internal 
feedback loop modeling is to consider ej, • • • , e, of E as new independent parameters and 
to find the worst possible bounds e, and e, for each e<; that is, e,- < e, < e, . This approach 
then reduces to the standard input-output decomposition problem with q independent 
parameters. However, e,’s may be functionally dependent to each other through actual 
parameter variations, Ap,’s. Whenever e,’s are closely related, this approach will result 
in a very conservative control design; furthermore, some valuable information on the 
structured parameter variations is not utilized in this approach. 

In order to exploite some structured or directional information on the plant parameter 
variations, the internal feedback loop gain matrix E is linearized about the nominal 
parameter set with respect to small Ap,-’s as follows: 

E = M\E x Ni (11a) 

A e = -MM 1 E 1 N 1 N (lib) 

where Ei contains only the actual, independent uncertain parameters. The standard 
form of an input-output decomposition such as Eq. (7), is obtained by re-defining M, 
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TV, and E as 


M «- MM, 

iV 4 — N\N (12) 

E *- E i 

In some cases, the plant parameter perturbations Ap,’s may possess a certain direc- 
tional relationship with the following form: 

A Pi = grf for t = 1 , * - • , / (13) 

where p, represents the direction and magnitude of A p,- and S is a scalar variable which 
represents the system uncertainty. Parameter variations involving a single parameter 
variable are referred to as uni- directional parameter variations, while those involving 
more than one parameter variable are referred to as multi-directional parameter varia- 
tions. 

The multi-directional parameter variations are characterized as: 


A pi = gijSj for t = 1, • • • , / and j = 1, • • • , r (14) 

where represents the direction and magnitude of A p,- caused by 6 : for multi-directional 
perturbations, and r is the number of independent uncertain parameters. 

For such cases with multi-directional parameter variations, the internal feedback loop 
gain E in Eq. (8) becomes a function of 6/ s: 


E = 


A e l (6) 
0 


0 

A e g (6) _ 


(15) 


where 6 = [<5 l5 • ■ • , 8 r ) T . Since E is nonlinear in Sj's in general, the linearized input- 
output decomposition can be applied here, as in Eqs. (11) and (12), to be incorporated 
in the robust control synthesis. 


Stability /Performance Robustness 

A robust Hoc full-state feedback control synthesis technique presented in this section 
exploits the internal feedback loop modeling concept and the Hoo control theory. This 
new robust control design methodology is summarized in terms of three theorems. De- 
tailed proofs of these theorems can be found in [8-10]. Development and application of 
robust H 0 o compensator synthesis can be found in [13]. 
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The parameter uncertainty model given by Eq. (7) and the nominal plant described 
by Eq. (2) can be combined as 


z p 


Gn Gn G\ 3 



z 

= 

G21 G22 G23 


W 

X 


1 

S 

0 

s 

co 

0 
1 


U 


w p = - Ez p 
u — — K x 


(16a) 

(16b) 

( 16 c) 


where w p and z p are, respectively, the fictitious input and output, E is the fictitious 
internal loop gain matrix, and K is a full-state gain matrix to be determined. 

The closed-loop system, but with the fictitious internal loop open, becomes: 



w p 



(17a) 

(17b) 


where 



T u = G u - G 1 Z K(I + G 33 K)~ 1 G 31 (18b) 

T n = G u - G 13 K(I + GmK)~'G 32 ( 18 c) 

T21 = G21 — G23JC ( X + <J 33 -^) - 1 G" 3 i ( 18 d) 

T 22 = G22 — G23K (I + G33K) * G 3 2 ( 18 e) 


The actual closed-loop transfer function matrix from w to z with plant perturbations 
becomes 


T, w = T 22 - T 21 E(I+ T u E)- 1 T 12 (19) 

Note that, in Eqs. (16) and (17), the parameter uncertainty does not appear in the 
transfer function matrices. Equations (17) can be used for the stability /performance 
robustness characterization. Sufficient conditions for robust stability and performance 
are provided by the following Theorems 1 and 2. 


Theorem 1 (Stability Robustness) 

T zw (s,aE ) Va € [0,1] is robustly stable for ||2?|| < e, and c > 0, if 

|[T„(s)|| 00 < e -1 

8 




where e is a measure of the magnitude of the plant parameter uncertainty E in Eq. (8). 


It is seen that T u determines the stability robustness with respect to parameter 
uncertainty. Small ||Tn||oo allows large parameter variations for closed-loop stability. 
For this reason T u is often referred to as robustness function [12]. The above theorem 
provides a sufficient condition for the closed-loop stability, resulting in a conservative 
control design. Since the condition in Theorem 1 is concerned with a deterministic 
bound, the control theory can be employed for the internal feedback loop model. 
The next theorem provides a sufficient condition for guaranteed performance robustness. 


Theorem 2 (Performance Robustness) 

T zw (s,aE ) Vq € [0,1] is stable, and ||T, u ,(s,aJE?)|| 00 < 7 Va € [0,1] with ||E|| < 7 -1 , 
if 

IITIU < 7 (20) 

where T and T zw are defined in Eqs. (17) and (19), and 7 is an upper bound for the 
desired performance specification. 


The above two theorems provide conditions for robust stability and performance of 
the perturbed closed-loop system in terms of Tn and T in Eqs. (17) and (18). The 
following re-definition of z, w , and the associated matrices enables us to employ the 
standard state-space representation given by Eq. (1): 


z <— 



w 


B x 


M x Bi 


C \ 


D n 


0 N w 
0 D n ’ 


D 12 


w p 

w ’ 
N x ' 
c i r 
N u ' 

D 12 


( 21 ) 


The following theorem [8] gives a robust H ^ controller which satisfies the condition in 
Eq. (20). 


Theorem 3 (Hoc Full-State Feedback Controller) Assume that 

(i) (A, B 2 ) is stabilizable and (C lt A) is detectable, 

(ii) Dl [ C, D„] = [0 I ], 
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(iii) the rank of Pn{ju) is m 2 for all u>, and 

(iv) D n = 0 . 

Given the above assumptions (i) through (iv), there exists an internally stabilizing 
controller such that, for the closed-loop transfer matrix T in Eqs. (17) and for a given 
design variable 7 , 

mice < 7 

if and only if the following Riccati equations 

0 = A t X + XA- X{B 2 Bl - ^B t Bj)X + C T X C X ( 22 ) 

have unique symmetric positive semi-definite solution X such that 

A - (. B 2 B\ - ^B x B t x )X and A - B 2 B*X are stable. 

A state-feedback gain that satisfies < 7 , where 7 is a design variable spec- 

ifying an upper bound of the perturbed closed-loop performance T zw , is then obtained 
as 

K = B\X (23) 

In order to achieve the desired closed-loop performance over all frequencies, T zw is 
often formulated to include frequency-dependent weighting matrices. (A proper selec- 
tion of the weighting matrices is an important step in any optimization-based design 
techniques, such as the linear-quadratic-gaussian (LQG) control and i/oo-optimization.) 
In this paper, constant diagonal weighting matrices are used. Inverses of the diagonal el- 
ements of the weighting matrix is referred to as weigthing factors. The weighting factors 
and 7 represent relative input-output levels and overall closed-loop performance level, 
respectively. In the Appendix, the usage of constant weightings, scaling, and orthogonal 
transformations on u, w, and z for practiced implementation of Theorem 3 are briefly 
summarized. 


3. Space Station Model 

The robust control synthesis technique developed in Section 2 is applied to the Space 
Station subject to large payload operations which cause significant changes in the mo- 
ments of inertia of the system. Dynamical equations of the Space Station are briefly 
reviewed (for details, see [3,4]). 
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The Space Station in a circular orbit is expected to maintain local-vertical and local- 
horizontal (LVLH) orientation during normal mode operation. For small attitude devi- 
ations from LVLH orientation, the linearized equations of motion can be written as: 


Space Station Dynamics: 


= n 


Attitude Kinematics: 


CMG Momentum: 


u>i ' 
tl>2 

U>3 

/31, 2/32, I 33 — 122 

—I 32 , 0, In 

I'll ~ All ~ 2/l2 > —Il3 


— ti8 3 = u>i 

&2 — n = U >2 

83 + nd\ = u>3 

hi — nh 3 = 
h 2 = u 2 

/13 + nh\ = U3 


u>i 

U^2 

UJ3 

'81 ■ 

82 

(24) 


(25a) 

(25b) 

(25c) 

(26a) 

(26b) 

(26c) 


■ /ll 

In 

1 13 


hi 

1 22 

1 23 


hi 

1 32 

I 33 . 



+ 3rr 


+ n 


I 33 ~ I 22) 
In, 

-I 13, 


Al, 

■f 33 — hi, 
—I 23, 


0 

0 

0 


5 

CM 

1 

1 


' -Ui + Wi 

CO 

CO 

+ 

— U 2 + U>2 

L -/« 


— U 3 + W 3 _ 


where (1, 2, 3) are the roll, pitch, and yaw control axes whose origin is fixed at the 
mass center, with the roll axis in the flight direction, the pitch axis perpendicular to 
the orbit plane, and the yaw axis toward the Earth; (0 1} 8 2 , #3) are the roll, pitch, yaw 
Euler angles of the body axes with respect to LVLH axes which rotate with the orbited 
angular velocity, n; (wj, w 2 , u> 3 ) are the body-axis components of the absolute angular 
velocity of the station; (In, I 22 , h 13 ) are the principal moments of inertia; I XJ ( i ^ j) 
are the products of inertia; (hi, /i 2 , /13) are the body-axis components of the CMG 
momentum; (ui, u 2 , U3) are the body-axis components of the control torque caused by 
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CMG momentum change; (w i, w 2 , u> 3 ) are the body-axis components of the external 
disturbance torque; and n is the orbital rate of 0.0011 rad/sec. 

Note that the products of inertia cause three-axis coupling as well as a bias torque in 
each axis. Fortunately, most practical situations with small products of inertia permit 
further simplification in such a way that pitch motion is uncoupled from roll/yaw motion. 
For the case where the control axes are nearly aligned with the principal axes (I x = In, 
h = In, and I 3 = 733), Eqs. (24) become 

thj -I- nkiu> 3 + 3n 2 ki0i = — 6i«x + biw x (27a) 

u> 2 + 3n 2 fc 2 0 2 — —bjU 2 + b?W 2 (27b) 

d> 3 - nk 3 ui = -&3U3 + b 3 w 3 (27c) 

WnPrP 

h={h-h)!h, 6i=i//i, 

* 2 =(/i-/ 3 )// 2 , b 2 =l/h, 

h=(l2-h)/h, h= l// 3 . 

Inertia matrices of the Phase 1 Space Station as well as the assembly flight #3 are 
listed in Table 1. In this paper, only the Phase 1 configuration is considered. The 
uncontrolled Space Station with such inertia properties is in an unstable equilibrium 
when 6i = 0 (i = 1,2,3). Also included are expected aerodynamic disturbances which 
are modeled as bias plus cyclic terms in the body-fixed control axes: 

w(t) = Bias + A n sin(nt + <f> n ) + A 2n sin(2 nt + <^ 2 n) (28) 

The cyclic component at orbital rate is due to the effect of Earth’s diurnal bulge, while 
the cyclic torque at twice the orbital rate is caused by the rotating solar panels. The 
magnitudes and phases of aerodynamic torque in each axis are unknown for control 
design. 

4. Space Station Control 

A robust JToo control design for the Space Station is described here. The Space Sta- 
tion is desired to have a control system which accommodates the periodic disturbances 
and large inertia variations. In [3,4], a periodic-disturbance accommodating controller 
is developed for the Space Station, and the disturbance rejection filters for the control 
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of hi, 62 , O 3 are assumed to have the following forms: 


&i + ( n ) 2 ai = h x 

(29a) 

A + (2n) 2 A = hi 

(29b) 

<22 H“ ( n ) 2 <>2 = A 

(29c) 

A + (2n) 2 /? 2 — A 

(29d) 

^3 + ( n ) 2q 3 = A 

(29e) 

03 + (2n) 2 A = A 

(29f) 


The pitch control logic, involving the single control input u 2 and eight states, is then 
expressed as 


U 2 = K 22*2 


(30) 


where K 22 is a 1x8 gain matrix and x 2 is the state vector defined as 


x 2 — [02 02 h 2 Pi 2 q 2 ®2 &2 Pip • 


(31) 


The CMG momentum and its integral are included to prevent CMG momentum build- 
up. 

Similarly the roll/yaw control logic is given by two control inputs, Ui and 113, and 
sixteen states: 

r l r Er' - -. fc'.. l r i 

(32) 




Ku 

*13 


*1 

. “ 3 . 


. *31 

* 33 . 


*3 


where K ^' s are 1x8 gain matrices and 


*1 = [$i hi Jhi q?i di fii fti] T , 


*3 = [03 w 3 ^3 P 1 3 a 3 d 3 /?3 fc ] 7 • 


(33a) 

(33b) 


Directional inertia variations for the Space Station are modeled as 


[a/, 

A7 2 

M3 

II 

►4? 

<d 

0 7, ] 

(34a) 

[ AA 

A7 2 

II 

<1 

7 2 7 3 ] 

(34b) 

[ A Ii 

A7 2 

■< 

•*0 

II 

< 

0 —h ] 

(34c) 

[ A h 

A7 2 

< 

II 

<1 

— 7 2 O] 

(34d) 

[ aa 

A7 2 

II 

•4? 

<1 

I 2 -h] 

(34e) 
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where <5,’s represent the amounts of directional parameter variations with respect to the 
nominal inertias I\, I 2 , and I 3 . The directional variation involving 6, is called a 5, -inertia 
variation in this paper. As discussed in [4], there exist physical bounds for Si’s due to 
the inherent physical properties of the gravity-gradient stabilization and the moments 
of inertia itself. Table 2 summarizes such physical limitations on Si -inertia variations. 
As discussed in [3,4], the Phase 1 Space Station becomes unstable for as little as — 7 % 
variation in I 3 and +8 % variation in Ji, because of the inherent physical nature of the 
problem. 

The robust controller synthesis in this paper is primarily concerned with the S 1- and 
^-inertia variations. In particular, the ^-inertia variation is physically caused by the 
translational motion of the payload along the pitch axis, as illustrated in Fig. 1. 


Pitch Control 

The pitch-axis dynamics with nominal inertias are described as: 

d_ 
dt 


'02 ' 


0 1 ' 

' e 2 ' 

^2 


— 3 n 2 k 2 0 

>2. 


+ 


0 
b 2 


[-“2] 


(35) 


where the external disturbance is not included since it is accommodated by the distur- 
bance rejection filter. 

Since the ^-inertia variation does not affect the pitch dynamics (i.e., k 2 and b 2 remain 
constant), only the ^-inertia variation, where only 62 has uncertainty, is considered for 
the pitch axis. 

An input-output decomposition of the perturbed control distribution matrix A B 2 in 
Eq. (35) is obtained as 


A B 2 = 




1 

h J 


= -MEN 


and 


M = 


0 

b 2 


E = S 2 , N = -1, 


where b 2 = l/I 2 for the nominal inertia and S 2 = 1/(1 + S 2 ) — 1. 

The fictitious input w p and the fictitious output z p for the pitch axis with the (52- 
inertia variation are then expressed as 


z p — N [— u 2 ] = «2 


(36a) 
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(36b) 


Wp — Zp Zp 


These equations replace the parameter variations in Eq. (35) as follows: 


d_ 

dt 


62 


*— H 

0 

1 

' 82 ‘ 

. ** . 


— 3 n 2 k 2 0 J 

. h . 


+ 
= u 2 


0 

in 


Wp - u 2 ] 


w„ 


— S 2 z p 


(37a) 

(37b) 

(37c) 


where k 2 = (h — h)/h with the nominal inertias. 

The robust control problem for Eqs. (35) now becomes a disturbance attenuation 
problem for Eqs. (37), to which Theorems 2 and 3 can be applied. Note that z p contains 
only the control input u 2 . This uncertainty in the control loop introduces a necessary 
tradeoff between stability robustness and performance. 

Equations (37) are now augmented by the pitch CMG momentum dynamics described 
by Eq. (26b) and disturbance rejection filters described by Eqs. (29c) and (29d). The 
augmented state vector is x 2 as defined in Eq. (31), and the controlled output z is also 
formed as 

z = [ * 2 

«2 J 

With proper selections of the weighting factors, scaling, and orthogonal transforma- 
tions, as discussed in the appendix, the augmented system equations are transformed 
to satisfy the assumption (ii) in Theorem 3. The performance specification bound 7 is 
chosen to be 1, and a set of weighting factors used in this paper is summarized in Table 
7. 

By solving the Riccati equation, Eq. (22), a robust Hoo full-state feedback controller 
for the pitch axis is obtained with a control gain matrix listed in Table 3. The closed- 
loop eigenvalues of the nominal system with this gain matrix are listed in Table 4. 
Stability margins of this new robust Hoo controller with respect to the inertia variations 
are compared in Table 5 to those of the previous LQR design in [3]. A significant margin 
of 70% for the <5 2 -inertia variation is achieved (compared to the 34% margin of the LQR 
design). As can be seen in Table 4, however, this new pitch controller has a closed-loop 
pole at — 8.29n which is relatively large compared to that of the conventional LQR design 
of Ref. 3. As discussed in [14], an Hoo controller often achieves the desired robustness 
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by having a high bandwidth for a single input system. The pitch axis design here is 
such a case; but the robust H control design for the multi-input case to be discussed 
in the next section has a remarkable stability robustness margin with nearly the same 
bandwidth as the conventional LQR design. 

Roll/Yaw Control 

Consider the roll/yaw dynamics with nominal parameters described by 

0 

■Til 

1 


' 9 X ' 

r 

Ai 


03 


(^3 

L 


+ 


0 


1 


n 

- 3 n 2 ^i 

0 

0 ■ 

— n 


0 

0 

0 

nk 3 

0 

■ 0 

0 ' 




61 

0 



-«i 

0 

0 



-Ms . 

. 0 

63 . 






’ 0i ' 


U>i 


03 


. ^3 . 


(38) 


where the external disturbances are not included. Variations in k x , k 3 , 6j, and 63 caused 
by perturbations in the moments of inertia are approximated as follows: 

ll 1 1 1 1 


A k, a -^ + M 

h h 1 z 


A 6, & - b j 


A h 


A63 — — 63 


Ah 


J 1 -> "•» J 

-*i h 

where k\ = (h — h)/I\, k 3 = (/ 2 — h)/h, h — l//i, and 63 = I//3, for the nominal 
inertias. 

In particular, for the < 5 i-inertia variation, the above parameter variations become: 

Ak 1 — — (&i -f* l)( 5 j, A b\ — — ^i^i) 

Ak 3 at - (t, + l)(x ar -63^. 

•*3 -*3 

An input-output decomposition of the perturbed system matrix A c is then obtained 


as: 


A e = -ME [ N x N u ] 
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or 




0 

0 ' 






- ■ 

M = 

bi 

0 


E 



0 < 



0 

0 




0 *1 




0 

b* . 







N x = 

f — 3n 2 /„ 


0 

0 

— n/ tt 




0 


0 

0 

5 


N u = 

-1 

0 1 





— 

0 

— 


1 





where /„ = 1? — I 3 4- 1\ and /(, 

<311 

— 1 \ + h- 





The fictitious input w p and output z p for the roll/yaw control design incorporating 
the inertia variation are expressed as 


z p = N x x + N u 


-Ui 

-it 3 


w p = -Ez p = -<$1 z p 


(39a) 

(39b) 


where x = [#x u>\ 63 u> 3 ] T . The perturbed system is then expressed by the nominal 
system and the internal feedback loop as 

x 
z p 

w p 

where u = [u x u 3 ] r , B = M and 

0 1 n 0 

— 3n 2 ki 0 0 — nki 

-n 0 0 1 

0 nk 3 0 0 



= Ax + B[w p — u] 

(40a) 

= N x x — N u u 

(40b) 

P. 

1 

II 

(40c) 


Similarly to the pitch-axis design, the standard state-space representation given by 
Eq. (1) can be constructed by redefining w and z. A roll/yaw gain matrix of the robust 
Hoo controller is listed in Table 3 for the particular weighting factors chosen as in Table 
7. The closed-loop eigenvalues of the nominal system with this robust Hoo controller 
axe listed in Table 4. Stability margins of this new controller are compared to those 
of other previous designs in Table 6. Similarly to the pitch control design, the robust 
Hoo controller for the coupled roll/yaw axes has significant improvement in stability 
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margins over the standard LQR design (e.g., the 73% margin over the 44% margin for 
the ^-inertia variation). Contrary to the pitch case with a single control input, however, 
the robust H<x, controller for the roll/yaw axes with two control inputs has a relatively 
low bandwidth! In fact, the roll/yaw closed-loop poles shown in Table 4 are very 
comparable to those of LQR designs in [3-5]. 

5. Discussions 

Major results and contributions of this paper are summarized in this section. A 
robust control synthesis technique presented in Section 2, which is primarily based on 
the results in [10,13] and the state-space formulation of the Hoo control theory in [8,9], 
further exploits the concept of linearized, directional variations of nonlinear, structured 
uncertain parameters. Applications of this approach to the full-state feedback control 
design problem of the Space Station with uncertain inertia property have resulted in the 
following interesting results: (1) For the pitch control with a single input, the stability 
robustness improvement with respect to the overall inertia increases has been achieved 
mainly by having a relatively high bandwidth controller and (2) The robust control 
design for the roll/yaw axis with two control inputs has achieved significant stability 
robustness over the LQR design, even with relatively low bandwidth. In other words, the 
concept of linearized directional parameter variation, combined with the standard 
control theory, has been shown to be a practical way for designing parameter-insensitive 
controllers. 

For roll/yaw control, the ^-inertia variation was considered in robust control 
design to accommodate the moments-of- inertia variations caused by the translational 
motion of a large payload along the pitch axis (see Fig. 1). Since the ^-inertia variation 
does not affect the pitch dynamics, the ^-inertia variation was considered for the pitch 
control design. It is also emphasized that the closed-loop system with this new robust 
Hoo controller is stable for ±73% (Sx-inertia variation and for ±70% 6 2 -inertia variation, 
compared to the ±44% and ±34% S 2 stability margins of a typical LQR design. 

6. Conclusions 

A robust control synthesis technique for uncertain dynamical systems subject to 
nonlinear, structured parameter perturbations has been presented, which is based on 
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the Hoc control theory and the internal feedback loop modeling concept. This technique 
was applied to the multivariable, full-state feedback control design problem of the Space 
Station, resulting in remarkable stability margins with respect to the moments-of-inertia 
uncertainty over the conventional linear-quadratic-regulator designs. The linearized, 
directional parameter variation concept was shown to be a proper way of accommodating 
the nonlinear, structured parameter variations in the design of a parameter-insensitive 
controller. 


Appendix 

In general, the Hoc control theory considers frequency-dependent weighting matrices 
for the shaping of closed-loop transfer function T IW . Proper selection of the weighting 
matrices, however, is not always obvious. One practical way is to use a constant diagonal 
weighting matrix and a normalized output equation. Proper scaling and orthogonal 
transformations can be employed to satisfy the assumptions in Theorem 3. 

Consider a system given by 


x = Ax + BiW + JB 2 « 

r> (2) u 

L 1 -'i2 J 

where D ^ is assumed nonsingular. 

Define r 2 ( i), r 2 ( 2 ), and r w be the weighting factors with dimensions of p 1} m 2 , and 
m l5 respectively. The weighting matrices Q, R, and W are then defined as: 

Q= [diag{r^<i) }] — 1 
R= [diagfryi)}] -1 
W = [diagfr^}] -1 



Define normalized variables as 


2 ^) = 
z (2) = 


w = W w 
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A scaling factor S for u is also defined as: 


u = Su 


Substituting the above new variables into the system equation gives 


= Ax + B\W 1 w + B-iS l u 


5 -1 u 


' Q-»SW ’ 


’ cS J) ' 


n (1) 1 

i/l2 

R-'zW 


0 

X + 

£j(2) 

xy 12 J 


The controlled output equation can be rewritten as 


r jtw 

- z (2) 


QC\ 

o 


a) 


x + 


with a QR decomposition of the matrix 

QD[\ ] 


rd { 3 


= P 


QD 

RD 


0 

L 


(1) 1 
12 

( 2 ) 

12 


s -1 « 


where P is an orthogonal transformation matrix and L is lower-triangular (or generally 
nonsingular). 

If the control scaling matrix S can be defined as 

S = X" 1 

the following system equation then satisfies the assumptions (ii) and (iv) in Theorem 3: 

x = Ax + J9 a W -1 «? + BiLu 


zW 

zW 


= P 1 


QC\ 

0 


(i) 


* + 


0 

I 


u 


where P T does not affect the H norm property. 

The system matrices axe redefined, to be implemented in a computer software (e.g., 
CTRL-C), as 

B X +-B X W- 1 , B 2 +-B 2 L, 


C, 


QC\ 

0 


(i) 


b 2 

D 12 


0 

I 


Finally the actual control gain matrix K is obtained by re-scaling the normalized gain 
matrix K as: 

K = S~ l K 
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The weighting factors selected for the example design of this paper are listed in Table 7 . 
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Table 1: Space Station Parameters 


Parameters 

Assembly 

Flight#3 

Phase 1 

Inertia (slug-ft 2 ) 


In 

23.22E6 

50.28E6 

1 22 

1.30E6 

10.80E6 


23.23E6 

58.57E6 

hi 

— 0.023E6 

-0.39E6 

hz 

0.477E6 

0.16E6 

1 23 

-0.011E6 

0.16E6 

Aerodynamic 

torque (ft-lb) for Phase 1 

w x 

1 + sin (nt) + 0.5sin(2nt) 

w 2 

4 + 2 sin (nt) + 0.5 sin(2 nt) 

U> 3 

1 + sin (nt) + 0.5sin(2 nt) 


Table 2: Physical bounds for ^-inertia variations 


Variation Type 

Lower Bound 

Upper Bound 


-78.5 %* 

00 

*2 

1 

*— • 
O 
O 

O 

$ 

00 

^3 

-2.3 %* 

+7.6 %» 

6 4 

-64.6 %* 

+16.4 % f 

* 

-2.1 %* 

+7.6 %* 


Mue to pitch open-loop characteristic. 

*due to roll/yaw open-loop characteristic. 

*due to triangle inequalities for the moments of inertia. 
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Table 3: Robust H «, controller gains for the Phase 1 Space Station 


Pitch JiTjj 

Roll/Yaw K t 

Units 


6.885E+2 

2.559E+2 

ft-lb/rad 


4.092E+5 

1.448E+5 

ft-lb-sec/rad 


2.648E— 3 

1.495E-3 

ft-lb/ft-lb-sec 


— 5.563E— 7 

4.948E— 7 

ft-lb/ft-lb-sec 2 


— 1.147E— 10 

— 4.142E-10 

ft-lb-rad 2 /ft-lb-sec 3 


5.193E-7 

2.263E— 7 

ft-lb-rad 2 /ft-lb-sec 2 


-9.374E-10 

-9.864E-10 

ft-lb-rad 2 /ft-lb-sec 3 


— 3.783E— 7 

-3.204E— 7 

ft-lb-rad 2 /ft-lb-sec 2 

4.531E+2 

1.800E+2 

4.115E+2 

ft-lb/rad 

2.607E+5 

8.914E+4 

3.719E+5 

ft-lb-sec/rad 

1.169E— 2 

5.124E-4 

2.015E— 3 

ft-lb/ft-lb-sec 

4.518E— 6 

— 3.149E-7 

— 1.997E— 7 

ft-lb/ft-lb-sec 2 

5.673E— 5 

— 1.567E— 5 

— 8.042E— 5 

ft- lb- rad /sec 2 

3.598E— 2 

— 6.513E— 2 

— 2.127E-3 

ft- lb- rad /sec 

— 1.722E— 5 

1.892E— 4 

— 2.489E-4 

ft-lb- rad/sec 2 

6.626E— 2 

8.709E— 4 

2.506E— 2 

ft-lb-rad/sec 


Table 4: Closed-loop eigenvalues of the Phase 1 Space Station with robust H <*, controller, 
in units of orbited rate, n = 0.0011 rad/sec 



Momentum/ 

Attitude 

Disturbance 

Filter 

Pitch 

-0.54 ± 0.54/ 

-0.10 ±1.05; 


-1.53, -8.29 

-0.10 ±2.03; 

Roll/Yaw 

-0.20, -0.21 

—0.13 ± 1.01; 


-0.31 ± 0.87; 

-0.33 ± 1.1 8; 


-0.82 ± 0.85; 

-0.10 ±1.99; 


-2.31 ±0.65; 

-0.27 ± 2.06; 
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Table 5: Pitch-axis stability robustness comparison 



LQR 

Robust Hoo 

% 

6 

6 

6 

6 

Si 

-99 

OO 

-99 

OO 

s 2 

-89 

34 

-99 

70 

S 3 

-17 


-27 

7 

s< 

-19 



16 

6s 

-30 


-31 

7 


'6 and 6 are lower and upper bounds, respectively. 


Table 6: Roll/yaw stability robustness comparison 


| 

LQR 

Local* 

Robust Hqo 

% 

6 

6 

6 

S 

S 

6 

Si 

-78 

44 

-64 

29 

-78 

73 

s 2 

-99 

43 

-67 

30 

-99 

71 

^3 

-61 

80 

-60 

61 

-58 

77 

Si 

-64 

64 

-64 

35 

-64 

99 

6s 

-51 

68 

-48 

50 

-49 

66 


* local feedback control (decentralized control). 

'6 and S are lower and upper bounds, respectively. 
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Table 7: Weighting factors used in the example design 


Weightings 

Pitch 

Roll 

Yaw 

Units 

e 

1.5 

1.7E-3 

2.1E-3 

rad 

6/u> 

3.8E-3 

6. IE— 7 


rad /sec 

h 

8.7E+3 

1.7E+4 


ft- lb- sec 

fh 

3.7E+4 

1.2E+5 

1.8E+5 

ft-lb-sec 2 

a 

3.5E+4 

1.6E+8* 

9.2E+2 

sec 2 /rad 

a 

2.1 

1.7E+5* 

8.0E-1 

sec/rad 

0 

2.5E+4 

3.7E+r 

1.5E+3 

sec 2 /rad 

0 

1.0 

3.7E+5* 

3.5E-1 

sec/rad 

u 

2.7E-2 

1.0E-1 

1.2E-1 

ft-lb 

Z P 

2.7E-2 

5.0E-2 

5.0E-2 

ft-lb 

w p 

1.0E-2 

1.0E-2 

1.0E-2 

ft-lb 


*in units of ft-lb-sec 3 /rad 2 , ft-lb-sec 2 /rad 2 , ft-lb- 
sec 3 /rad 2 , and ft-lb-sec 2 /rad 2 , respectively. 
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